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Abstract. We report on the results of the first XMM-Newton satellite observation of the luminous and helium-rich O-type 
subdwarf BD-r37° 1977 carried out in April 2014. X-ray emission is detected with a flux of about 4x10“^^ erg cm“^ s“^ 
(0.2-1.5 keV), corresponding to a /x//boi ratio 10~^; the source spectrum is very soft, and is well fit by the sum of two 
plasma components at different temperatures. Both characteristics are in agreement with what is observed in the main-sequence 
early-type stars, where the observed X-ray emission is due to turbulence and shocks in the stellar wind. A smaller but still 
significant stellar wind has been observed also in BD-l-37° 1977; therefore, we suggest that also in this case the detected X-ray 
flux has the same origin. 
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1. Introduction 


Among the hot subdwarf (s d) stars, whic h are evolved He- 
core burning low-mass stars ( Hebei 2009h . the sdO stars are 
those which show the highest temperatures (Teff > 40 kK). 
Apart from this characteristic, sdO stars (sdOs) are character¬ 
ized by a wide range of values for the surface gravity (log(g) 
= 4-6.5) and helium abundance (-3.5 ^ log(nHetT-H^) ~ 3)' 
Indeed they form a rather heterogeneus class of stars, which in¬ 
cludes both He-po o r and He-rich stars ( Heber & JeffervlfT992 ; 
Heber et al. 2006 : Hirsch et~^ 2008 ). and ‘luminous’ and 


‘compact’ stars, according to their low or high surface grav¬ 
ity, respectively (lNat)iwotzkill2008l) . This variety of properties 
is probably the consequen ce of different evolutionary histories 
(lHebeill2009l: lGeieill20 1 -5h : in the case of the compact stars, the 
He-poor ones are post-EHB stars, while the origin of the He- 
rich ones might be either the merging of two He-core or C /O- 
core white dwarfs ( Iben 1990l : Saio & JeffervI 200C , 2002 ) or 
the so-called late hot-flasher scenario ( Brown et alJ 2001 ): in¬ 
stead the luminous sdO stars are post-AGB stars. Evolutionary 
models suggest that most sdO stars are the outcome of the evo¬ 
lution of single stars, but some of them could descend from 
binary systems which underwent a common-envelope phase; 
in this case it is possible that the sdO stars has a compact com¬ 
panion, typically a white dwarf (WD). 


Up to now, sdO stars have been deeply investigated in the 
optical/U V domain, where several of them are rather bright; on 
the other hand, only few of them are known as X-ray sources. 
In the case of the binary HD 49798, the detection of pulsed 


(P = 13.18 s) soft X-rays (llsrael et al.lll997h indicates that 
this emission originates from accr etion onto a compact ob ject, 
most likely a massive white dwarf ( Mereghetti et al.ll2009l) . Eor 
this binary we detected an evident X-ray emission also when 
the compact companion is eclipsed by the sdO star, suggest- 
ing the possibility of in trinsic X-ray emission of the sdO star 
(IMereghetti et al.ll2013l) . Another sdO star recently detected at 
X-rays is BD h- 37° 442: the XMM-Newton observation of this 
He-rich star revealed soft X-ray emission, with a spectrum sim¬ 
ilar to that of HD 49798, and a possible periodicity of 19.16 s 
(at 3 cr confidence level), wh ich suggests that also BD -1- 37° 442 
has a compact companion (ILa Palombara et al.N2012h . In or¬ 
der to enlarge the sample of sdO stars observed at X-rays, we 
performed with Chandra HRC-I a survey of a complete flux- 
limited sample of sd O stars and discovered th ree additional 
X-ray emitting stars ( La Palombara et ^ 20141): the luminous 
and He-rich sdO star BD-i-37° 1977 djefferv & Hamannll2010h 
and the compa ct (logft/) > 6) and He-poor stars Eeige 34 an d 
BDh- 28° 4211 dTheill et alJll99ll:IZanin & Weinbergeill 19971) . 


In this paper we report on the results of a follow-up obser¬ 
vation of BD -1-37° 1977, performed with XMM-Newton, which 
allowed us to investigate in detail the spectral and timing prop¬ 
erties of the X-ray emission discov ered with Chandra . This 
star was identified as an sdO star by IWolff et al.l (119741) . who 
detected several emission He lines but no H lines in its blue 
spectrum. Their spectroscopic analysis gave a surface gravity 
logp ^ 4.5 and a temperature T ^ 50 kK; comparable values 
(T ~ 55 kK and logp ~ 4.0) were estimated from the low res- 
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Table 1. Main parameters of the sdO stars BD +37° 1977, BD +37° 442, and HD 49798. 


Parameter 

Symbol 

BD+37° 1977 
Value Reference 

BD +37° 

Value 

442 

Reference 

HD 49798 

Value Reference 

Surface gravity 

logs 

~4.0 

1 

4.00 ± 0.25 

4 

4.35 

6 

Luminosity (Lq) 

L 

25,000 

1 

25,000 

1 

14,000 

6 

Effective Temperature (K) 

TeS 

48,000 

2 

48,000 

2 

46,500 

6 

Magnitudes 

U 

8.67 

3 

8.57 

5 

6.76 

7 


B 

9.93 

3 

9.73 

5 

8.02 

7 


V 

10.17 

3 

10.01 

5 

8.29 

7 

Distance (kpc) 

d 

~2.7 

2 


4 

0.65 ± 0.10 

8 

Terminal wind velocity (km s“^) 

Voo 

2,000 

2 

2,000 

2 

1,350 

9 

Mass-loss rate (M© yr“^) 

M 

Iq-8.2 

2 

10-8.5 

2 

10-8.5 

6 


References: 1 - [Darius e t alj 1 19 791): 2 -l^fferv & H amann H2010l); 3 _ _ _ _ _ 

iHamannl boiol) : 7 - Landolt & Uomotol J2007 ): 8 - Kudritzki & Simon ( 1978h : 9 - iHamann et al.l 1 198lh 


Jordi et al.l h99lh: 4 -iBauer & Husfell ( Il995h : 5 - iLandold l ll973l) : 6 ■ 
iHamani 


olution lUE spectrum dParius et al.lll979h . which also gave an 
estimate of the star luminosity (log(Lboi/ Lq = 4.4). There is no 
evidence for a compact companion for BD +37° 1977. The pos- 
sible detection of an infrared excess (at 2-tT confidence level, 
Ulla & Theilll ( 1998h l. if confirmed, could imply a companion 
star of spectral type earlier than G4. Based on high-resolution 
ultraviolet and optical spectra and on UV-optical-IR photome¬ 
try, and using the latest gener ation of models for spheric ally ex¬ 
panding stellar atmospheres, Jeffery & Hamann ( 2010l) found 
that BD+37° 1977 is characterized by a significant mass-loss 
rate (M = Mq y“^). Their analysis yielded revised and 

better constrained spectral parameters, very similar to those of 
BD +37° 442, the other X-ray detected He-rich sdO star. The 
best-fit distance modulus obtained from this analysis is DM = 
12.2, corresponding to a distance of ~ 2.7 kpc. Our observa¬ 
tion of BD+37° 1977 with Chandra HRC-I provided a detec¬ 
tion with a count rate CR = 3.6 ^q g cts s“^, which, assuming 
a spectrum similar to that of HD 49798 and BD +37° 442, im¬ 
plies an X-ray luminosity Lx ~ 10^^ erg s“^. For comparison, 
in Table [T] we report the main parameters of the three luminous 
sdO stars detected in X-rays. 


2. Observations and data analysis 

BD +37° 1977 was observed with XMM-Newton in April 2014. 
At that time the source was visible only for the first 20 ks of 
each XMM-Newton orbit; therefore, five different observations 
were performed, between April 14 and 24 (see Tabl e |2l). The 
three EPIC cameras, i.e. one pn ( Striider et alJl200l1) and two 
MOS ( Turner et~^ 200 ih . were always operated in full frame 
mode, with time resolution of 73 ms for the pn and of 2.6 s for 
the two MOS cameras; taking into account all the observations, 
the total effective exposure time was, respectively, of ~ 34.5 ks 
and ~ 50 ks. For all cameras the medium thickness filter was 
used. 

We used version 13.5 of the XMM-Newton Science 
Analysis System (SAS) to process the event files. For the data 
analysis we selected only the events with pattern in the range 
0-4 (i.e. mono- and bi-pixel events) for the pn camera and 
0-12 (i.e. from 1 to 4 pixel events) for the two MOS. For 
each camera, we merged together the data of the five obser¬ 


vations and accumulated the images in various energy ranges. 
We found that BD +37° 1977 is significantly detected at the co¬ 
ordinates R.A. = 9^ 24”^ 26.4", Dec. = +36° 42' 52.8", which 
differ by 0.7" from the position of BD +37° 1977. This differ¬ 
ence is consistent with the ^ 2" r.m.s. astrometric accuracy of 
XMM-Newtot^ In each of the five observations a point source 
is clearly detected below 0.5 keV, while considering the five 
observation merged together the source is detected up to ~ 1.5 
keV (Fig.[I]i. This implies that the source spectrum is very soft. 
All the observations were partly affected by high instrumental 
background. However, since the spectrum of the instrumental 
background is rather hard, the background contamination has 
a limited impact on the source spectral analysis: therefore, we 
considered the whole data set, without rejecting the time inter¬ 
vals with the highest particle background. The source net count 
rate in the 0.15-1.5 keV range is (1.7 ± 0.2)x 10“^ cts s“^ and 
(2.5 ± 0.3) X10“^ cts s“^ for the pn and each of the two MOS, 
respectively. 

For the timing and spectral analysis, we used the data of 
the whole observation and the three EPIC cameras; we ex¬ 
tracted the source events from a circular region with radius 
15" centered at the source position, while the corresponding 
background events were accumulated on circular areas free of 
sources and radii of 30" and 120" for the pn and the two MOS 
cameras, respectively. We converted the arrival times to the 
solar system barycenter, then we combined the three datasets 
in a single event list. The background-subtracted light curve 
of BD +37° 1977 does not show any variability on time scales 
from hundreds of seconds to the observation length. We looked 
for possible periodicities in the X-ray emission, but we found 
no evidence of periodic signals; this search was unsuccessful 
not only for the five individual observations, but also when con¬ 
sidering together all of them. In all cases we estimated an upper 
limit of ^ 30 % on the pulsed fraction, for a sinusoidal modu¬ 
lation between 1 and 5000 s. 

For the spectral analysis we considered first the pn data, 
since the source soft spectrum and the lower sensitivity of the 
MOS cameras at low energies strongly reduced the count statis¬ 
tics. We verified that the addition of the MOS data gave consis- 

’ http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018.ps.gz 
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Table 2. List of observations of BD +37° 1977 performed by XMM-Newton. 


Revolution 

Observation ID 

Start time 

(YYYY-MM-DD hh:mm:ss) 

Effective exposure 
pn (ks) MOS (ks) 

2627 

0740140301 

2014-04-14® 15:00:25 

4.2 

8.5 

2629 

0740140501 

2014-04-18@14:43:21 

7.0 

11.3 

2630 

0740140401 

2014-04-20® 15:45:37 

7.9 

10.0 

2631 

0740140601 

2014-04-22® 14:36:30 

9.9 

14.4 

2632 

0740140701 

2014-04-24® 14:06:44 

5.5 

6.8 



Fig. 1. EPIC pn image of the sky region around BD +37° 1977 in the energies ranges 0.15-0.5 {left), 0.5-1.5 (centre), and 1.5-10 
keV (right). The green circle (15" radius) indicates the source position. 


tent results. We generated the response matrix and ancillary file 
using the SAS tasks rmfgen and arfgen. To ensure the ap¬ 
plicability of the statistics, the spectrum was rebinned with 
a minimum of 30 net counts per bin; then, we fitted them us¬ 
ing XSPEC (V 12.7.0). We only used the energy range 0.2-1.5 
keV since at higher energies the background dominates and the 
source flux is negligible. In the following, all the spectral un¬ 
certainties and upper limits are given at the 90 % confidence 
level for one intere sting parameter, and we assume a source 
distance of 2.7 kpc (Jeffery & Hamann 20101) ; we adopted the 
results of Anders & Grevessel ( 19891) for the solar abundances 
of the atomic elements. 


The source spectrum is very soft and we tried to describe 
it with different models (see Tabled. The fit with an absorbed 
power law (PL) is formally acceptable (x^ < 2) but gives a 
very large and unrealistic photon index (L ~ 5), while a fit 
with a blackbody model is rejected by the data (x^ > 2). An 
absorbed power law plus blackbody gives a good fit, but with 
unrealistic (for the power-law photon index) or unconstrained 
(for the blackbody normalization) values of the model parame¬ 
ters. We note that while this model is physically motivated for 
BD +37° 442, where the observed X-ray flux can be attributed 
to accretion onto a compact companion, this is not the case 
for BD+37° 1977, for which no evidence of a compact com¬ 
panion has been found. For this reason, we consider in the 
following the possibility that the X-ray emission detected in 


BD +37° 1977 has the same origin of that observed in the nor¬ 
mal, giant and supergiant early-type O stars. 


For a large sample of this type of stars observed with 
XMM-Newton the spectrum can be described by the sum of 
different thermal plasma components (MEKAL in XSPEC) , 
with temperatures between ~ 0.1 and ~ 5 keV (lNazal2009l) . 
Therefore, we tried to use the same approach also in the case 
of BD +37° 1977. We clearly found that, assuming solar abun¬ 
dances, with this model it was not possible to obtain an accept¬ 
able fit, even if we considered the sum of two MEKAL compo¬ 
nents at different temperatures (Fig.|2]and TableO. Therefore, 
we modified the model abundances by taking into account 
the values of the single chemical elements considered by 


Jeffery & Hamann ( 2010l) . Since there are no abundance mea¬ 


surements for BD +37° 1977, they adopted th e same overabun¬ 
dance of He, C, N, Si, and Fe obtained by IBauer & Husfeld 
(119951) for BD +37° 442, which is very similar from the spec¬ 
troscopic point of view. 


Assuming these abundances, a single thermal component 
provides an acceptable fit (xf < 2) but leaves large residuals at 
the high energies. Hence we considered a model composed by 
the sum of two absorbed components, at different temperatures. 
We checked that the estimated interstellar absorption is negli¬ 
gible and consistent with 0. Therefore we fixed it at A^h = 10^° 
cm“^, which is the total absorption value across the Galaxy in 
the direction of BD+37° 1977. In this way we found a good 
fit (xl = 0.71 for 13 degrees of freedom. Fig. |3ll with kTi = 
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0.2 0.5 1 2 


Energy (keV) 

Fig. 2. Top panel: EPIC pn spectrum of BD +37° 1977 with the 
best-fit model composed by the sum of two thermal plasma 
components, with solar abundances. Bottom panel: residuals 
(in units of a) between data and model. 



0.2 0.5 1 2 

Energy (keV) 


Fig. 3. Top panel: EPIC pn spectrum of BD +37° 1977 with the 
best-fit model composed by the s um of two thermal plasma 
components, with abundances from Llefferv & HamannI (I 2 OIOI 1 . 
Bottom panel: residuals (in units of cr) between data and model. 


Table 3. Summary of the best-fit parameters of BD +37° 1977 
obtained with different spectral models. 


Parameter 

Unit 

Value 

Power law 


iVn 

r 

xl 

cm ^ 

(3.5+2 7)x 10^° 

5.0^11 

1.76 

Degrees of freedom 

- 

14 

Blackbody 
A^’h cm“^ 

kT eV 

(0.9+^-^) xlO^® 
67t^^ 

xl 

- 

2.46 

Degrees of freedom 

- 

14 

Power law + Blackbody 
iVn cm-2 (2.0j:^'°)xl02i 

Tpl - 

kTsB 

eV 

Z,Z_i2 

xl 

- 

0.76 

Degrees of freedom 

- 

12 

Mekal + Mekal 


(with solar abundances) 

Nh 

cm“^ 

(1.8±1.2)xl0^° 

kTi 

kT2 

eV 

eV 


xl 

- 

2.15 

Degrees of freedom 

- 

12 

Mekal + Mekal 


Iwith abundances from Jeffierv&Hamann ('201011 

Na 

cm“^ 

1X10^° (fixed) 

kTi 

eV 

120±30 

kT2 

eV 

8401®®“ 

xl 

- 

0.71 

Degrees of freedom 

- 

13 


Mekal + Mekal 


(with abundances from Jeffier^&jlamMin (201^) 
(with free Oxygen abundance) 


Na 

cm“^ 

1X10^° (fixed) 

kTi 

eV 

lOOtfo 

kT2 

eV 


Oxygen Abundance 

- 

310±250 

xl 

- 

0.64 

Degrees of freedom 

- 

12 


120±30 eV and kT 2 = eV. The total flux in the en¬ 

ergy range 0.2-1.5 keV is /x = (d.Oj^g'g) xl0“^^ erg cm“^ 
s“^; although 78 % of the flux is due to the low-temperature 
component, the high-temperature component is significant at 3 
cr confidence level. The measured flux corresponds to a source 
luminosity Lx = (S.Slgg) xlO^^ erg s“^. The fit leaves some 
residuals at ^ 650 eV, thus suggesting the presence of the 
O VIII emission line (Pig. |2l. This could be due to an under¬ 
estimate of the real Oxygen abundance, since in our model we 
fixed it at the Solar value. Therefore we repeated the spectral 
fit with the same model but leaving the Oxygen abundance free 
to vary. In this way we obtained an improvement of the spec¬ 
tral fit (Pig. nil and we found that the best-fit abundance value 
is 310±250 times the Solar value; although it is not well con¬ 


strained, this is consistent with the expected Oxygen overabun¬ 
dance in BD+37° 1977. 


3. Discussion 


The XMM-Newton observation of BD+37° 1977 enabled us 
to constrain the flux and spec trum of the X-ray emissio n re¬ 
cently discovered by Chandra ( La Palombara et ^12014 ). The 
measured flux /x — 4 x 10“^^ erg cm“^ s“^ confirms the 
estimate provided by the Chandra detection. It implies a lu¬ 
minosity Lx — 3.3 X 10^^ erg s“^. Since the b olometric lu¬ 
minosi ty of BD+37° 1977 is Lh ni — 25,000 Lq ( Darius et al 
1979[ Llefferv & Hamanri 201ot) . the corresponding ratio is 
Ly^lLhoi ^ 10“®'®. This value is consistent with the ‘canoni- 
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Fig. 4. Top panel: EPIC pn spectrum of BD +37° 1977 with the 
best-fit model composed by the su m of two thermal plasma 
components, with abundances from Jeffery & HamanrJ ( 2010ll 
and free Oxygen abundance. Bottom panel: residuals (in units 
of a) between data and model. 


cal’ relation Lx 10“^ x Lboi obtained for the normal, giant 
and supergiant early-ty pe O stars, which are known since a long 
ago as X-ray source s ( Pallayicini et al. 198lt Sciortino et al 


199(i Ir.iidel & Na7dl2009l) .' The hypothesis of intrinsic origin 


for the X-ray emission of BD+37° 1977 is further supported 
by the spectral analysis. In fact, considering the likely possibil¬ 
ity of non-Solar composition, the spectrum can be successfully 
described by the sum of two thermal plasma components, as in 
normal early-type stars. 

It is interesting to compare the properties of the three sdOs 
for which X-ray spectral information is ayailable (Table S. 
The spectrum of HD 49798 during the eclipse phase can be 
described by the sum of three thermal plasma components 
(IMereghetti et al.ll2013h . While the hottest component is re¬ 
quired to account for the significant emission aboye ^ 2 
keV, the temperatures of the two coldest components are yery 
similar to those of BD+37° 1977. The 0.2-10 keV luminos¬ 
ity of HD 49798 during the eclipse phase is Lx ~ 3 x 
10^° erg s“^, i.e. one order of magnitude lower than that of 
BD +37° 1977. Howeyer, due to its lower bolometric luminos¬ 
ity, also for HD 49798 the X-ray/bolometric luminosity ratio 
is Lx/Lboi ^ 10“^. The spectrum of the other X-ray emit¬ 
ting sdO star, BD +37° 442, can also be fit with a similar the- 
mal model with the He and m etal abundances deriyed for this 
star djeffery & Hamann 2010 ) and the temperatures indicated 
in Table|4] Its flux /x — 6.4 x 10“^"*^ erg cm“^ s“^ implies a 
luminosity Lx — 2.9 x 10^^ erg s“^ (for a distance of 2 kpc), 
hence a X-ray/bolometric luminosity ratio Lx/Lboi = 10“®'^. 

These results indicate that the X-ray emission from these 
three luminous sdO stars is similar to that of normal O-type 
stars, which haye luminosities up to a few 10®® erg s“^. In these 
stars the X-ray emission is due to turbulence and shocks in their 
strong embedded winds ( Lucy & White 1980l : Owocki et al 


19881). Since the mass-loss rate in the radiation-driyen winds 


of early type stars scales with the bolometric luminosity, and 


the X-ray emission originates in the stellar wind , a correlation 
betwe en Lx and Lboi is not surprising (see, e.g.. lOwocki et al 


(l2013h ). In this respect it is interesting to note that the three X- 
ray emitting sdOs are among the few hot s ubdwarfs for which 
eyide nce of mass loss has been reported (IJeffery & Hamann 
2010h . Our results indicate that, eyen if the winds of sdO stars 
are rather weak (e.g. M = 1 0~®'^M(n yr~^ for BD+37 ° 1977 
according to the estimate of Jeffery & Hamann ( 2010h L they 
can produce X-ray emitting shocks, as in more luminous O type 
stars. In this framework, we note that our findings for sdO stars 
are supported by the m ethodology used and the results obtained 
by ICohen et alJ (120141) . who inyestigated the X-ray spectra of 
O-type stars with yery low mass-loss rates. 

The latest Chandra detections of sdO stars reinforce the hy¬ 
pothesis that also this type of stars is a class of X-ray sources; 
in this respect, the detection not only of the luminous stars 
but also of the compact ones is yery promising. In addition, 
we also note that in all cases the estimated Lx/Lboi ratio 
agrees with that found in the heayier early-type stars: there¬ 
fore also in the compact sdO stars the X-ray emission could 
be attributed to turbulence and shocks in their winds. In or¬ 
der to proyide an oyeryiew of all the sdO stars obseryed at 
X-rays up to now, in Fig. |5] we report, as a function of their 
bolometric magnitude, the X-ray flux of the detected stars 
and its upper limit for the undetected ones; for HD 49798, 
BD +37° 442, and BD +37° 1977 the flux yalue is based on the 
spectral fit proyided by the XMM-Newton data, while for the 
other sources it is based on the count rate yalue or limit pro¬ 
yided by Chandra HRC-I, assuming an emission spectrum sim¬ 
ilar to that of BD +37° 1977. For comparison, the two dashed 
lines trace the region corresponding to the typical X-ray-to- 
bolometric flux ratio for the normal O-type stars. The plot 
shows that almost all the obseryed stars are within this region, 
therefore the possible presence of intrinsic X-ray emission also 
from the stars undetected so far cannot be excluded. 
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Table 4. Summary of the best-fit parameters of the three sdO stars observed with XMM-Newton, when their spectrum is described 
with multi-temperature thermal-plasma components (MEKAL in XSPEC). 


Source 

kTi 

(keV) 

kT2 

(keV) 

kTa 

(keV) 

log(Lx/Lboi) 

Reference 

HD 49798 (in eclipse) 

0.13±0.02 

0 7] 

^•'^-0.19 

5 (fix) 

-7.1 

Mereshetti et al. 120131 

BD +37° 442 

0 I'7+0.02 
'-0.01 

fv '^'^+0.22 
'^-0.10 

- 

-6.3 

This work 

BD+37° 1977 

0.13±0.01 

0.79±0.15 

- 

-6.5 

This work 



Bolometric magnitude 


Fig. 5. Level of the X-ray flux (or its upper limit for the un¬ 
detected sources) of the sdO stars observed at X-rays, as a 
function of their bolometric magnitude. The upper and lower 
blue lines (corresponding to fx/fhoi = 10“®-^ and /x//boi = 
10“^'^, respectively) include the range of expected values for 
the main-sequence early-type stars; the red line correspo nds to 
/x//h ni = 10“®’^, which is the best-fit relation found by Nazd 
( 2009h for this type of stars. 
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